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Several puzzles about the data at high pr in heavy-ion collisions are 
listed. The resolution of them all is given in the framework of parton 
recombination. More specifically, it is the recombination of the soft and 
semi-hard shower partons that enhances the region 3 < pr < 9 GeV/c, 
and gives rise to the large p/ir ratio in AuAu collisions. The Cronin effect 
can be explained in terms of final-state interaction for both it and p. The 
structure of jets produced in AuAu is different from that in pp collisions. 
The suppression of Rc p in forward production can also be understood by 
extending the same hadronization scheme at r\ — to r\ > without the 
introduction of any new physics. 

PACS numbers: 25.75.-q 

In the past two years several features of the high-p^ data obtained by 
various experiments at RHIC are puzzling, and may be regarded as anoma- 
lies according to the "standard model." By standard model I mean that 
which has been standard in the treatment of hadron production at high pr, 
namely: a hard scattering of partons, followed by a fragmentation process 
that leads to the detected hadron. What I plan to show in this talk are 
evidences that all those anomalies can be resolved when the process of hard 
parton fragmentation is replaced by the recombination of soft and shower 
partons. The basic reason why the fragmentation model has worked so 
well for high-pT processes in leptonic and hadronic collisions, but poorly for 
heavy-ion collisions, is that there is a large body of soft partons in the latter 
case, but absent in the former. The hadronization of those soft partons by 
recombination with the semi-hard partons results in a significant enhance- 
ment in the intermediate-pr region that is missing in the fragmentation 
model. 
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If hard partons fragment in vacuum, whether or not they have lost energy 
while in transit in the dense medium, the fragmentation products should 
be independent of the medium. Thus the ratio of produced hadrons, when 
all else are the same, should depend only on the ratio of the fragmentation 
functions (FF), D(z). Given a parton, whether a quark or a gluon, its 
FF for the production of a proton D p (z) is much smaller than that for a 
pion D 7T (z). The observed data reveal several anomalies according to that 
picture. 

Anomaly 1. The ratio of proton to pion, R p / % , in Au+Au collisions is 
approximately 1 at pt ~ 3 GeV/c. 

Anomaly 2. The nuclear modification factor, Rcp, in d+Au collisions is 
greater for p than for ir at px ~ 3 GeV/c. 

Anomaly 3. The jet structure, i.e., the distribution of particles associated 
with a trigger, is different for jets produced in Au+Au collisions compared 
to that for p+p collisions. 

Anomaly 4- The azimuthal anisotropy parameter V2 is larger for baryons 
than for mesons for pt > 2 GeV/c. 

Another irregularity. Forward-backward asymmetry. 

Time and space do not permit all the topics above to be addressed 
adequately. I give only a sketch here. 

How can recombination solve the puzzles? First of all, let it be under- 
stood on general grounds that when a multi-parton state is to hadronize, it 
is far more efficient for a q and q to recombine than for a higher momentum 
q to fragment, assuming that the parton distribution is falling rapidly in 
momentum. That is simply because recombination involves the addition 
of two lower momenta of q and q, where the densities are higher, whereas 
fragmentation involves first the creation of a parton at higher momentum 
(at a cost in yield), and then the production of a hadron at some momen- 
tum fraction at the cost of another factor of suppression. The comparison is 
meaningful only when there are many soft partons moving collinearly with 
a hard parton, which is the case for heavy-ion collisions, but not for leptonic 
and hadronic collisions. 

The fragmentation process makes use of a phenomenological FF because 
it describes the non-perturbative process of hadronization that cannot be 
calculated in pQCD. Thus D(z) represents a black box, in which there are 
gluon radiation, quark pair creation, etc., that generate a shower of partons 
before hadronization. Although the density of shower partons cannot be 
calculated from first principles, those partons are nevertheless there, and 
their momentum distributions can possibly be determined phenomenologi- 
cally. That is what we have done in the framework of recombination pP, in 
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which we write 

xD(x) = [ — — F q q>(xi,x 2 )RM(xi,X2,x), (1) 

J X\ X2 

where F q qi is the joint distribution of a shower quark q and a shower an- 
tiquark q' that recombine to form a meson M. Rm is the corresponding 
recombination function. The FF on the LHS is known from phenomenolog- 
ical analysis of leptonic and hadronic processes. Rm is known from previous 
work on the recombination model. So F qq -i can be determined. If we assume 

that F q qi has the factorizable form SfSf , where i labels the hard parton 
that initiates the shower, then there are 5 types of Sj to be determined from 
5 types of Di, where i takes on the species: u, d, s, g, and j is allowed to be 
u,d,s, but not g, on the grounds that gluon conversion to qq relieves the bur- 
den of considering direct hadronization of gluons. The distributions Sj (x) 
of the shower partons have been determined in at a fixed Q = 10 GeV/c. 
The Q 2 evolution of Sj (x) was not considered, although it constitutes an 
interesting project in its own right. On the basis that hadron production in 
the intermediate pt region at RHIC depends crucially on the recombination 
of soft and shower partons, but not sensitively on the virtuality of Sj (x) we 
have proceeded to the study of the consequences of considering the shower 
partons in heavy-ion collisions 2 , and found some remarkable results. 
For pion production at large pt the inclusive distribution is 

dN w f dpi dp2 

p 



f JPiJElp ' q -,(p 1 ,p 2 )R n (p 1 ,p2,p), (2) 
dp J pi p 2 



where 



R?r(Pl,P2,P) = +P2 ~P)- (3) 



Similar equations can be written for proton production [2j. The essence 
of recombination is then in what one should include for F u j in case of 7r + , 
say, and for F UU( i for p. If we denote thermal parton distribution by T and 
shower parton distributions by <S, then they ought to be 

F u<I = rr + rs + ss, (4) 

F uud = TTT + TTS + TSS + SSS, (5) 



where the pure T terms give the soft component, and the pure S terms 
recover the fragmentation component. It is the mixed terms involving both 
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T and S that are new and dominate the intermediate-pr region, as we shall 
show below. 

To proceed, we need to specify T and S. For T we do not rely on any 
low-pT model, but determine it phenomenologically from the low-p? data 
of pion production, using the TT term of @ in ©. In that way we can 
attribute the enhancement at higher px directly to the TS term without 
raising any question on the reliability of the \ow-px model. For S we shall 
use the distributions Sj determined in 1 1 , convoluted with the distribution 
of hard parton i in Au+Au collisions, and then sum over i. More specifically, 
we write 

T(pi) = C Pl e-^ T , (6) 



5fe) = eE / dkkh{k)si(p 2 /k) , (7) 

i 

where fi(k) is the hard-parton distribution. C and T are parameters to be 
varied to fit dN^/pxdpT for px < 2 GeV/c. £ is the average fraction of hard 
partons that can get out of the dense medium to produce showers. Without 
discussing the details that can be found in [2], let me just show the results. 
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Fig. 2. Comparison of calculated p/ir ratio with data from 



Fig. 1 shows the pion distribution that exhibits the dominance of the 
TS component in the px region between 3 and 9 GeV/c. The sum of 
all components agrees well with data [3]. Similarly, the TTS and TSS 
components dominate over other components of the proton spectrum in 
the same pr region. The p/ir ratio is shown in Fig. 2; indeed, it reaches 
the level of 1 at around px ~ 3 GeV/c, as observed |3]. Without the 
thermal-shower recombination, the proton spectrum would be too low in the 
intermediate-^ region, so the first anomaly in the fragmentation picture 
is now satisfactorily resolved. It should be mentioned that the large p/ir 
ratio has also be obtained in other recombination/coalescence models using 
slightly different schemes to implement the hadronization process El • 

The second anomaly concerns the Cronin effect in d+Au collisions. The 
traditional explanation of the effect is that multiple scattering of partons in 
the initial state leads to the broadening of the px distribution of hadrons 
at high px in the final state. If those hadrons are produced by the frag- 
mentation of high-pr partons, then Rqp for protons should be much lower 
than that for pions. However, the data at RHIC reveal R^p > Rqp ^ or 
1 < px < 3 GeV/c [7j. This anomaly can be well explained in the recom- 
bination model when the consideration given above for Au+Au collisions is 
extended to d+Au collisions [B1IH1, as shown in Fig. 3. Although the soft 
partons in d+Au are not thermal in the sense of Au+Au, the soft com- 
ponent nevertheless plays a similar role and the recombination of soft and 
shower partons gives rise to components that can reproduce the tt and p 
spectra in px without any initial parton broadening. Thus it is the final- 
state rather than the initial-state interaction that is mainly responsible for 
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Fig. 3. Comparison of calculated ratios for Rcp for it and p with data from 0. 



the Cronin effect. The reason is simply that the formation of proton by the 
recombination of 3 quarks, some soft some semi-hard, is far more effective 
than the fragmentation of a hard parton. 

The third anomaly according to the fragmentation picture is that the jet 
structure in Au+Au collisions is different from that in p+p collisions. Data 
from STAR show that the charge multiplicity and total scalar pt in the 
near-side jets are significantly higher for Au+Au than for p+p, even when 
the trigger (that is the same for both) is included ^Uj. If the trigger were 
not included, one would expect the associated particle distributions to be 
drastically different for the two cases. Such behaviors cannot follow from the 
fragmentation of hard partons, once the trigger momentum is fixed to be the 
same. This problem has been studied in the framework of the recombination 
model, in which at least two shower partons in a jet must be considered, 
one for the trigger, the other for the associated particle. Specifically, for 
7r + (trigger) and tt + (associated) in a jet the 4-parton distribution has the 
structure 



777r" r 7T~ r 

^4 



(TS)(TS) + (TS)(SS) + (SS)(TS) 



(8) 



where the first pair of parentheses in each term correspond to the trigger, 
the second pair the associated particle. We have omitted the term (SS) (SS) 
in that equation because it is negligible in Au+Au collisions; however, it is 
the only term that is important in pp collisions. This point clearly reveals 
the difference between jets produced in heavy-ion and hadronic collisions. 
The difference becomes even greater when other types of associated parti- 
cles are included, since the thermal environment in heavy- ion collisions helps 
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the formation of other mesons and baryons in conjunction with shower par- 
tons. In Fig. 4 we show the distribution associated with a ir + trigger when 
7r + ,7r~ and p in the jets are all included The data ^U] are for all 

charged hadrons in both the trigger and the associated particles, and are 
therefore not exactly what we have calculated. Nevertheless, the agreement 
is remarkably good. 
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Fig. 4. Transverse momentum distribution of tt + ,it and p associated with a tt + 
trigger in central Au+Au collisions. Data are from JU| for all charged hadrons. 



The fourth anomaly concerns elliptic flow where vi for baryon exceeds 
that for meson. This phenomenon has nicely been explained by the coa- 
lescence model |12| . and the scaling of v% with the number of constituent 
quarks remarkably verified by the STAR data |13| . 

The final issue to be mentioned here is about the production of interme- 
diate-^ hadrons at large forward rapidity in d+Au collisions. It is the 
region where high hopes have been raised for the verification of a signature 
of color glass condensate. As with pQCD, the hadronization mechanism 
is fragmentation. BRAHMS data already show that Rcp at rj = 3.2 rises 
no higher than 0.5 at px ~ 3 GeV/c 14 j. This suppression is regarded as 
evidence for gluon saturation |15| . However, before novel physics is invoked, 
it is reasonable to ask whether the phenomenon can be understood in the 
conventional way, i.e., by extrapolating what is known to work at midrapid- 
ity to the forward region. We have preliminary results that show a general 
agreement between the data and the expectation from parton recombination 
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at all rj and pr- The spectra at forward rapidities are suppressed because 
there are less soft partons as r] is increased, resulting in less hadrons formed 
that rely on soft partons recombining with shower partons. At rj = 3.2 there 
are so few hard partons that most hadrons are the result of soft-soft recom- 
bination. We have also studied the backward-forward asymmetry and found 
that there is no need for a transition of basic physics from multiple scatter- 
ing in the initial-state interaction on the r] < side to gluon saturation on 
the rj > side. 

Our emphasis on the hadronization process in the final state provides 
a universal framework for the description of particle production at all r\ 
and pt, at all centralities. In that framework of interpreting the existing 
high-pT data from RHIC there are no features that are puzzling. In a sense 
that may be disappointing, since exciting new physics usually comes with 
anomalies. However, it is far better to have no puzzles than to be misled 
by false anomalies. 

All the theoretical results reported here were done in collaboration with 
C. B. Yang. This work was supported, in part, by the U. S. Department of 
Energy under Grant No. DE-FG02-96ER40972. 
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